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ABSTRACT 


This  computer  program  computes  a  ballistic  trajectory  from  burn  out 
or  any  point  thereafter  to  impact.  It  assumes  point  mass  and  a  spherical 
earth.  The  program  takes  into  account  air  resistance,  the  variation  of  the 
gravitational  and  centrifugal  field  with  altitude,  and  earth's  rotation.  It 
uses  a  system  of  co-ordinates  (X,  Y,  and  Z)  rigidly  connected  with  the 
rotating  earth. 

The  computer  program  is  written  in  FORTRAN  II  and  an  autocoder 
deck  has  been  compiled.  The  program  is  run  on  the  IBM  1410  computer. 

It  is  being  adapted  for  use  on  the  IBM  1620  computer  for  use  by  the  Army. 
This  program  may  be  obtained  from  the  Scientific  Programming  Section; 
Program  No.  SP9-  It  was  programmed  by  W.  B.  Warren  of  the  Scientific 
Research  Staff. 
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SYMBOLS 


Symbol 


Computer 

Symbol 


p 

BETA 

Bearing  angle  -  Clockwise  horizontal  angle  from 
north  to  direction  of  flight 

6 

DELTA 

Flight  path  angle  -  Angle  between  missile  axis 
and  local  horizontal.  Angle  is  positive  when 
missile  is  climbing. 

<t> 

PHI 

Latitude  -  Positive  in  northern  hemisphere, 
negative  in  southern 

e 

THETA 

Longitude  -  Increases  to  the  east.  When  firing 
east,  0Q  -  0,  to  the  west,  0Q  =  ir/2 

X 

X 

See  Figure  1  -  Co- ordinate  system  with  center 
of  earth  as  origin.  Z  is  the  polar  axis.  X  and 

Y 

Y 

Y  are  in  the  plane  of  the  Equator  and  are  90° 
apart.  The  X  axis  corresponds  to  the  point 

Z 

Z 

where  0  ==  0,  All  are  positive  as  shown. 

Y 

GAMMA 

Range  angle  -  Vertex  at  center  of  earth  measured 
from  the  point  where  data  were  put  into  the  program 
to  the  position  attained  along  the  trajectory. 

Ax 

AX 

Acceleration  in  the  X  direction 

Ay 

AY 

Acceleration  in  the  Y  direction 

Az 

AZ 

Acceleration  in  the  Z  direction 

A 

t 

AT 

Acceleration  along  the  trajectory 

vx 

VX 

Velocity  in  the  X  direction 

VY 

VY 

Velocity  in  the  Y  direction 

Vz 

VZ 

Velocity  in  the  Z  direction 

Vt 

VT 

Velocity  along  the  trajectory 

Ra 

RA 

Distance  from  center  of. earth  to  re-entry  body 

m 


Symbols  (Cont. ) 

Computer 

Symbol  Symbol 


H  ALT 

R  R 


r 


Distance  from  surface  of  earth  to  re-entry  body 
Radius  of  earth  -  Constant  (20,  902,  890  ft.) 
Geocentric  location  vector 


So 

g 

U 


c r 

u 


GO  Acceleration  of  gravity  at  sea  level  -  constant 

(32  .174  ft/sec2) 

GRAV  Acceleration  of  gravity  at  altitude 

U  Earth's  rotation  in  radians  per  second  - 

constant  (  2tt/24  x  3600) 

SIGMA  Angle  of  a  spherical  triangle  -  see  Figure  3 

AMU  (gcR2)  ^ 

BALCO  Ballistic  coefficient  -  Estimated  from  shape 

and  dimensions 


p  DEN  Air  density  from  tape  for  IBM  1410  program  - 

computed  internally  in  IBM  1620  program 

PHYSICAL  CONSTANTS 


R  =  20,  902,  890  feet 
go  =  32,  174  ft/sec2 
U  =  2tt/24  x  3600  radians/sec. 
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INTRODUCTION 


This  computer  program  may  be  used  to  study  the  free  flight  and 
re-entry  portions  of  a  ballistic  trajectory.  The  information  calculated 
and  printed  by  the  computer  includes  time  of  flight,  location  in  both 
X,  Y,  and  Z  from  earth's  center  and  latitude,  longitude  and  altitude  from 
earth's  surface,  acceleration,  velocity,  flight  path  angle,  bearing  angle 
and  range  in  nautical  miles. 

The  model  chosen  for  use  in  this  computer  program  was  taken  from 
References  1  and  2.  Point  mass  and  a  spherical  earth  are  assumed.  The 
basic  acceleration  equations  are  developed  in  Reference  1.  Reference  2 
discusses  a  similar  computer  program  which  uses  the  same  acceleration 
equations  and  many  of  the  other  relationships  used  in  this  computer  program. 

A  portion  of  the  computer  program  discussed  in  Reference  3  was 
used  to  compute  density  in  the  IBM  1620  Computer  program.  This  portion 
is  an  integral  part  of  that  program. 

Library  Functions  Used: 

SORT 

SIN 

COS 

EXPF 

ATAN 
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ANALYSIS 


Initial  Conditions 

The  computer  program  will  accept  initial  conditions  at  burn  out 
or  at  any  time  thereafter  until  impact.  The  items  of  information  needed 
are  punched  on  two  IBM  cards  which  are  placed  at  the  end  of  the  deck. 


Items 

6 

P 

0 

<t> 

H 

R 


U 


Comments 

Positive  when  missile  is  climbing. 

Clockwise  horizontal  angle  from  north 

Initial  0Q  =  0  when  flying  eastward  and  0O  =  ir/2 
when  flying  toward  the  west 

Positive  in  northern  hemisphere 

Altitude  above  earth's  surface 

20,  902,  890  ft.  (earth's  radius) 

Velocity  along  the  trajectory 

32.174  ft/sec2 

2tt/24  x  3600  radians  per  second 
Ballistic  coefficient 


All  angles  are  put  into  the  computer  in  radians.  The  computer 
program  is  designed  to  accept  and  compute  English  units.  All  angles  are 
printed  in  degrees.  (Program  SP-9-E) 

The  computer  program  may  be  modified  to  accept,  compute  and  print 
metric  units  by  using  metric  units  for  H,  R,  Vt>  gQ  and  by  using  conversion 
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factors  for  the  altitude  and  density  in  the  density  subroutine  in  the  IBM  1620 
computer  program.  A  deck  in  metric  units  has  been  prepared.  (SP-9-M) 
Development  of  Equations 


In  order  to  provide  the  proper  initial  conditions,  the  following  set 
of  co-ordinate  transformations  were  used  as  an  integral  part  of  the 
computer  program.  (See  Figure  1): 


X  =  Ra  cos  <)>  cos  3  (1) 

'l 

Y  =  Ra  cos  <|>  sin©  (2) 

Z  =  Ra  sin  <f>  (3) 

By  definition: 

H  =  R  -  R  (4) 

a 

It  follows  that: 

Ra  =  (X2  +  Y2  +  Z2)^  (5) 


In  order  to  convert  the  input  velocity  along  the  trajectory  to 
components  in  the  X,  Y,  and  Z  directions,  the  angles  6,  (3,  4>,  and  0  are 
used.  See  Figure  1 .  6  is  the  angle  the  vehicle  makes  with  the  local  hori¬ 

zontal. 


Vt  sin  6  =  a  velocity  component  along  Ra  (6) 

Vt  sin  6  sin  <(>  =  a  velocity  component  along  the  Z  axis  (7) 

Vj.  sin  6  cos  <|>  =  a  velocity  component  along  Ra  projected  (8) 

into  the  plane  of  the  Equator 

Vt  sin  6  cos  <(>  sin  0  =  a  velocity  component  along  the  Y  axis  (9) 
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sin  6  cos  (|>  cos  0  =  a  velocity  component  along  the  X  axis  (10) 

Vj.  cos  6  -  a  velocity  component  parallel  to  the  earth's  surface  (11) 

V  cos  6  sin  p  =  a  velocity  component  along  a  parallel  of  (12) 

*  latitude 

Vt  cos  6  sin  (J  sin  0  =  a  velocity  component  along  the  X  axis  (13) 

Vt  cos  6  sin  p  cos  0  =  a  velocity  component  along  the  Y  axis  (14) 

Vj.  cos  6  cos  p  =  a  velocity  component  along  a  meridian  of  (15) 

longitude 

Vt  cos  6  cos  p  sin  (j>  =  a  velocity  component  along  Ra  projected  (16) 

into  the  plane  of  the  Equator 

■  Vt  cos  6  cos  p  sin  <|>  sin  0  =  a  velocity  component  along  the  (17) 

Y  axis 

V  cos  6  cos  p  sin  <(>  cos  0  =  a  velocity  component  along  the  (18) 

X  axis 

Vt  cos  6  cos  p  cos  <|>  =  a  velocity  component  along  the  Z  axis  (19) 

The  following  three  equations  are  integral  parts  of  the  computer 
program.  The  signs  (  +  or  -)  of  the  sin  and  cos  functions  of  6,  p,  tj> ,  and 
0  were  analyzed  as  the  angles  changed.  The  proper  sign  was  given  to  the 
components  used. 

Combining  (10),  (13),  and  (18): 

=  Vj.  (sin  6  cos  <j>  cos  0  -  cos  6  sin  p  sin  0  -  cos  6  cos  p  sin  <j>  cos  0)  (20) 

Combining  (9),  (14),  and  (17): 


Vy  =  Vt  (sin  6  cos  si*1  9  +  cos  6  sin  p  cos  0  -  cos  6  cos  p  sin  <J>  sin  0)  (21) 


Combining  (7)  and  (19); 


V"  =  V  (sin  6  sin  4>  +  cos  6  cos  (3  cos  §) 
z  t 


The  component  velocities  can  be  used  to  find  the  velocity  along  the 


trajectory. 


V  =  (V  2  +  V  z  +  V,2)2 
t  x  y  z 


If  one  takes  into  account  the  air  resistance,  the  variation  of  the 
gravitational  and  centrifugal  field  and  earth's  rotation,  the  equation  of 


motion  for  the  ballistic  vehicle  becomes: 


where 


=  mg  +  W  +  k  +  k 


mass  of  vehicle 


radius  vector  from  center  of  earth  to  vehicle 


acceleration  due  to  gravity 


W  =  air  resistance 


Coriolis  force 


centrifugal  force 


Examining  equation  (24)  term  by  term,  one  obtains: 


Ra  Ra 
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where 


R  =  radius  of  earth 

R  =  distance  from  earth's  center  to  vehicle 
a 

W  =  ■  vcv11  *  (26) 

where  c  and  n  are  factors  which  depend  on  air  density,  the  shape  of  the 
vehicle  and  the  velocity; 

kc  =  2m  (v  x  U)  (27) 

where 

U  =  angular  velocity  of  the  rotation  of  earth  in  radians  per  second 

and 

_  i 

kn  =  mU2  r  cos  cf>  =  mU2  (X2  +  Y2)2  (28) 


In  this  analysis,  the  air  resistance  was  assumed  to  be  proportional 
to  the  square  of  velocity  so  that  equation  (26)  becomes 


—  v  p  g  V  m 

W  =  — - 

2 


(c^) 


(29) 


Writing  equations  (25)  through  (29)  in  component  form,  we  obtain; 


u  X  pVtgVx 

A  =  -  JT  -  — /  — r  +  2  v  u  +  u2  x 
x  Ra  i  w  \  y 


(c^) 


(30) 


A  =  _  n  Y  pVtg  Vy 

y  R  J 


(-sU 

VcdA/ 


-  2  V„U  +  U2  Y 


(31) 
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(32) 


A„ 


u  Z 

Rj 


P  Vt  g  Vz 

2(^) 


Integrating  these  equations  in  computer  form,  we  obtain 


\ 

=  Vx(n  - 

*  1) 

+ 

A  At 

X 

(33) 

VYn 

=  VY(n  - 

■  1) 

+ 

Ay  At 

(34) 

% 

=  VZ(n- 

■  1) 

+ 

Az  At 

(35) 

Xn  = 

Xn-  1 

+ 

vx 

At 

(36) 

Y  = 
n 

Yn-1 

+ 

Y 

y 

At 

(37) 

Z  = 
n 

Zn-1 

+ 

V 

z 

At 

(38) 

From  Figure  1,  the  relationships  pertaining  to  the  angles  are: 


<(>  =  arc  sin 


0  =  arctan 


6  =  arc sin 


Ra 

Y 

X 

ARa/At 

vT 


where  ARa/At  =  average  rate  of  change 


t  in  altitude 

In  spherical  trigonometry,  the  law  of  cosines  states: 
cos  a  =  cos  b  cos  c  +  sin  b  sin  c  cos  A 


(39) 

(40) 

(41) 


(42) 
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The  law  of  sines  states: 

sin  A  _  sin  B  _  sin  C  1^2) 

sin  a  sin  b  sin  c 

where  angles  A,  B,  and  C  represent  the  angles  of  the  spherical  triangle  on 
the  earth's  surface  and  a,  b,  and  c  represent  the  sides  of  the  triangle  and 
the  angles  subtended  by  them  at  earth's  center.  See  Figure  2. 

Now  refer  to  Figure  3. 

cos  y  =  cos  Co  4>0  cos  Co<f>  +  sin  Co<t)0  sin  Co  4>  cos  A0  (44) 

where 

A0  =  0  -  0Q  and 

Co  <|)  =  90  0  -  >  then  cos  Co<j>Q  =  sin  (j>0  and  sin  Co4>0  =  cos  4>0 

Co<j>  =  90°  -  <j>  »  then  cos  Co<j>  =  sin  4>  and  sin  Co<j>  =  cos  <(> 

Equation  (44)  then  becomes 


cos  y  =  sin  <()0  sin  <(>  + 

cos  <j>o  cos  (j)  cos  A0 

(45) 

y  =  arccos  (sin  <t>0  sin  <j> 

+  COS  <j)Q  cos  cj>  cos  A  0  ) 

(46) 

Then  range  =  yR  ,  y  is  in  radians.  Earth's  radius,  R  =  3440.  239  nautical 
miles . 

Range  =  y  (3440.23  9) 

Refer  again  to  Figure  3. 

sin  cr  _  sin  A9 
sin  Co4>0  ~  sin  y 

cos  <|>n  sin  A0 

sin  <r  =  - - - 

sin  y 


(47) 


but  sin  Co4>o  =  cos  cj>0 


(48) 


(49) 


z 


FIGURE  3 


or  = 


(50) 


arc  sm 


-( 


COS  <j>. 


sm 


sin  A6 

V 


) 


Sigma  (<r)  is  used  to  compute  bearing  (P).  In  order  to  minimize  the 
error  caused  by  assuming  that  the  trajectory  follows  an  exact  great  circle 
(it  would  follow  an  exact  great  circle  if  the  earth  were  not  rotating),  a 
smaller  spherical  triangle  is  used  in  the  computer  program. 

Let  symbols  with  subscript  1  stand  for  the  next  to  last  values 
computed.  The  symbols  without  subscripts  stand  for  the  newest  values 
computed.  Let  y3  equal  the  range  increment  angle. 

Equation  (46)  now  becomes: 


y3  =  arccos 


sin  4>i  sin  <(>  +  cos  4>i  cos  <t)  cos  (0  -  0i) 


(511 


Equation  (50)  becomes 

.  /cos  <bi  sin  (6  -  0i )  \  /c->, 

<r  =  arcsinl- _ If _ : _ Li  )  (5 2) 

v  sin  y3  / 

The  computer  uses  the  tangent  function  to  determine  the  value  of 
sigma  (<r) .  The  value  determined  will  lie  between  --90°  and  +90°.  A  graph 
showing  sigma  plotted  against  beta  is  shown  in  Figure  4. 

From  Figure  4  we  can  see  that  if  beta  is  between  0  and  90° 

(3  =  <r  (53) 

If  beta  is  between  90°  and  270°, 

p  =  it  -  <r  (54) 
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If  beta  is  greater  than  270 
p  =  2tt  +  <r 


(55) 


The  five  equations  (51)  through  (55)  are  integral  parts  of  the 
computer  program. 

In  order  to  compute  At ,  which  is  the  acceleration  along  the 

trajectory,  let  V  =  the  last  velocity  computed  along  the  trajectory, 

Vt  j  =  the  next  to  last  velocity  computed  along  the 
trajectory 

At  =  the  time  increment  being  used  the  computer 

then 


At 


(56) 


DISCUSSION 
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p 

!L 


<— ikJ 
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The  computer  program  assumes  a  spherical  earth  and  a  point-mass. 
Other  conditions  are  those  which  actually  exist  at  earth's  surface  and  at 
altitude  from  the  surface. 

The  time  increment  used  by  the  computer  is  one  second  during  the 
portion  of  the  trajectory  whe re  conditions  change  slowly.  During  the 
re-entry  portion  of  the  trajectory  when  the  altitude  of  the  re-entry  body  is 
below  150,000  feet,  the  time  increment  is  decreased  to  0.  1  second  to 
increase  accuracy.  In  order  to  give  desired  accuracy,  this  program  was 
written  in  double  precision.  The  program  uses  about  39,  950  digits  of 
memory.  It  is  believed  that  this  program  produces  highly  accurate  results. 

For  some  of  the  angle  relationships,  it  is  assumed  that  the  trajectory 
follows  an  exact  great  circle.  The  errors  caused  by  this  assumption  are 
very  small  because  of  the  short  time  increments  used  and  the  methods 
chosen  to  compute  these  relationships.  These  relationships  are  computer 
outputs  and  are  not  used  to  drive  the  computer. 


wen 
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FLOW  CHART  FOR  FREE  FLIGHT  AND  RE-ENTRY  TRAJECTORY 
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Density  Constants 
for  Zone  1 


Density  Constants 
.  for  Zone  2 


ALT-40,  000 


|  Density  Constants 
for  Zone  4 


Density  Constants 
for  Zone  5 


ALT-175,  000 


ALT-270,  000 


Density  Constants 
for  Zone  6 


Density  Constants 
for  Zone  7 


RO  = 


t _ 

<R3D+  R2)(D+  R1)(D+  RO) 
e _ 

100, 000 


ALT  -  400,00 


RO 

DEN  =  GRAY 


r  > 

7  98 

V  J 
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FORTRAN  STATEMENTS 
English  Units 
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CARD  TO  PRINT  80/80 


FREE  FLIGHT  AND  RE-ENTRY 
PROGRAMMED  BY  W.B. WARREN 
FOR  RESEARCH  STAFF 
10  NOVEMBER  1962 
PROGRAM  NO.  SP  9  ENGLISH  UNITS. 


PHYSICAL  CONSTANTS 

GO  H  GRAVITATIONAL  ACCELERATION  AT  SEA  LEVEL  -  32.174  FT/SEC**2 
U  H  EARTH  ROTATION  -  X2*P I n/£24*3600n  RADIANS/SEC 
R  H  RADIUS  OF  EARTH  -  20,902,890  FEET 
BALCO  It  BALLISTIC  COEFFICIENT  -  CONSTANT 
MU  *  *  2  U  G0*R**2  H  3  2.174*20,902,890**2 

OTHER  INPUT  DATA 

DEL  //ANGLE  BETWEEN  HORIZONTAL  AND  MISSILE  AXIS — POSITIVE  UP 
BETA  //AZIMUTH  ANGLE  CLOCKWISE  FROM  NORTH  TO  MISSILE  AXIS 
THETA  //LONGITUDE  OF  MISSILE 
PHI  //LATITUDE  OF  MISSILE 

H  //HEIGHT  OF  MISSILE  ABOVE  EARTH  SURFACE 

RHO  //ATMOSPHERIC  DENSITY  --TABLE 
RA  H  MISSILE  ALTITUDE  FROM  CENTER  OF  EARTH 
VO  H  MISSILE  BURNOUT  VELOCITY  //VT  AT  T//0 
RE A  C  100, DEL, BETA, THETA, PHI , ALT 
REAC100, R, VO, GO, U, BALCO 
READ  100, T, CELT, DELT2 
CONO  //  0. 

SENSE  LIGHT  1 
CON  15  //  0. 

CONI  //  0. 

AACC  //  0. 

C0NTR//0 . 

AABB  II  1. 

BET  A 1 #  BETA  -  .001 
ALT9  #1000. 

RN//G. 

VT  II  VO 
VT  1  //VO 

AMU  //  SQRTFSG0*R»*2n 

IPAGE  U  1 

BET  AO  II  BETA 

THETO  H  THETA 

PHIC  H  PHI 

RA  H  ALTER 

OVX  //VO*SINF3DELn*COSF$PHIu*COSF$THETAa— VO*COSF%DELo»COSF!JBETAa*SIN 

lF%PHIu*COSF^THETAn-VO*COSF^DELa*SlNF%BETAD*SINFlTHETAa 

OVY  //VO*SINF?DELa*COSFXPHIa*SINFUHETAD-VO*COSF%DELn*COSF*BETAa*SIN 
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CARD  TO  PRINT  80/80 


1FSPH  Io*S  INFSTHETAn£VO*COSFS!DELa*S  INFSBE  T  An*COSFSTHET  An 
0 \IL  #VQ*S  tNFSDELn*SINFSPHIn£VO*COSFSnELQ*COSFSBETAn*CO$F:SPHIa 
X  #RA*COSFSPHIn*COSF*THETAu 
Y  #RA*COSFSPHI a»S INFSTHETAn 
Z  #  RA*SINFSPHIn 
RA1  «  ALT  C  R 
PRINT  4999, IPAGE 
PRINT  5000 
PRINT  5001 
PRINT  5027 
PRINT  5002 
PRINT  5003 
PRINT  5004 
PRINT  5010 
PRINT  5011 
PRINT  5012 
PRINT  5013 
PRINT  5014 

PRINT  5015, DEL , BETA , THETA , PHI 
PRINT  5016, ALT,BALC0,RA,V0 
GAMMA  H  0. 

2462  GRA  V  II  S 32 . 174 *R*« 2a/R A ** 2 
IF  SALT  -  ALT9a900,901 ,901 

900  IF  SCON  1 5  -  1 . n3840  »  901 ,901 
3840  PRINT  5050 

5050  FORMAT  % 1HK , 3 5 X , 6HZ EN I THn 
CON  15  //  1. 

901  IF  SAL T-400GOO.n4200, 4201, 4201 
4200  IF  SALT  -  4CC00.nl, 2, 3 

1  D  #  ALT/IOOCO. 

R3#  -0. 37728875 E- 02 
R  2  #  0.52352523E-02 

Rl#  -0.31047929E800 
R011  0.89444960E801 

GO  TO  50 

2  CONTINUE 

3  IF  SALT-80000. Q4, 5,6 

4  CONTINUE 

5  0  II  SALT-40000. o/lOOOO. 

R3#  0.2329860 4 E- 05 

R2H  0.21620268E-03 
R1/H0.47974207E&00 
ROH  0. 754024 19EC01 
GO  TO  50 

6  IF  SALT  -  160000. n7, 8, 9 

7  0  II  SALT-80000. n/10000. 

R3#  0.2635981  IE-04 

R2 #  0 . 8094  3560E-02 

Rl#  -0.51509227EC00 
RO#  0. 56272847EG01 
GO  TO  50 

8  CONTINUE 

9  IF  SALT-175000. nlO, ll,  12 
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10  CONTINUE 

11  0  //SALT- 160000.  d/IOOOO. 


R3W 

-0.57 

066660E-04 

R2W 

0.30080CCQE-03 

Rl# 

-0.36 

384 1 7  3£  800 

ROW 

0.20 

423758E601 

GO 

TO  50 

12 

IF 

SALT- 

270000. nl3, 

14, 

15 

13 

0  # 

ULT 

-175000. n/L 

0000. 

R3  # 

-0.8 

4014648E-03 

R2# 

-0.1 

6797257E-02 

Rl# 

-0.3 

2931657E800 

ROW 

O'.  1 

502 1266E801 

GO 

TO  50 

14 

CONTINUE 

15 

IF 

%  ALT- 

290000. al6, 

17, 

18 

16 

CON 

TINUE 

17 

0  W 

%AIT 

-27C000.n/l 

GOOO . 

R3# 

0. 

R2# 

0.36 

480C09E-03 

R 1 W 

-0.61 

376240E800 

ROW 

-0.25 

094864E801 

GO 

TO  50 

18 

IF 

SIALT- 

350000. nl9, 

20, 

21 

19 

0  U 

-290000. a/1 

0000. 

R3W 

0.1 

7070042E-02 

R2  # 

-0.  1 

65792 18E-02 

Rl# 

-0.6 

389  87 1 8E  800 

ROW 

-0.3 

7292327E801 

GO 

TG  50 

20 

CON 

TINUE 

21 

0  # 

%  AL  T 

-350000. a/l 

OOOO. 

R3W 

-0.3 

8550824E-02 

R2W 

0.6 

590 1 94  5E-01 

Rl# 

-0.6 

8705559E600 

ROW 

-0.7 

258956 3 E801 

50  R0W3EXPFmR3*DSR2n*D6R  lo*D8RQnn/ 1 00000. 

OEN  U  RO/GRAV 

79 80 AX #2. *VY*UCU**2*X-X%AMU**2*Xd/RA**3d-*XDEN*GRAV  »VX*VTo/£2. *BALCO 

1  DD 

OAV  WU**2*Y-2  .*VX*U-$3AMU**2*YD/RA**3ci-UDEN*GRAV  *VY*VTn/X2.*BALC0 
Inn 

AZ#l-l.n*!«%?:AMU**2*Z^/RA**3aG*X0EN*GRAV  *V Z*VTa/12 .  *BALCOaaa 
IF  %C0N8o4000, 699, 1569 
699  CON 8  H  C0N8  8  1. 

GO  TO  601 
A 20 1  OEN  W  0. 

GO  TO  798 

1569  IF  25ALT  - 1 50000 . a  14 34 ,  1 434 , 3 1 04 
1434  IF  %AACCn4000, 2460, 2461 

2460  A ACC  H  AACC  8  1. 

GO  TO  3104 

2461  CELT  #  0  E  L  T  2 
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3104  T  It  T  l  DELT 
422  THETl  H  THETA 
BETA2  It BETA1 
BET  A  1  ft  SETA 
SIGMA10SIGMA 
PH  II  It  PHI 

704  PHI  #  ATANFSZ/*SQRTF*RA**2-Z*»2oon 
THETA  //ATANFSY/Xo 
OELRA  It  *RA  -  RAla/OELT 

DEL  ft  ATANF**DELRAo/*SQRTF*VT *«2-DELRA**2nan 
421  OTHETtfThETA  -THETl 

GAMMA3#  ATANF«XSQRTF*l.-**SINF*PHll'a*SlNFXPHIaaC»COSF*PHI'lD*COSFXP 
lHIo«C0SF*DTHETann*#2do/*SINF*PHl ln*SINFSPHlacCOSF*PHlla*COSF*PHln* 

2CQSF*DTHET  non 

SIGMA  it  ATANF3i%SINFa!0THETQ*C0SF*PH[  lnn/*SQRTF«*SINF*0AMMA3aa*»2-*S 

lINF%OTHETn*COSF*PHI  inci#*2ana 
IF  *1.5688  -  *ABSF*SIGMAaaa4099, 4160, 4150 
4099  IF  SSIGMA1O4102, 4101, 4102 

4101  OELSIG  H  .001 
GO  TO  4120 

4102  DELSIG  II  'SIGMA  -SIGMA1 

4120  IF*BETAl-BETA2n4128, 4130, 4129 

4128  BETA#  BETA  1-* A BSF*D EL  S I  Goa 
GO  TO  4300 

4130  BETA#  8ETA1G.001 
GO  TO  4300 

4129  BETA#  BETA  1  C*ABSF*DELS IGaa 
GO  TO  4300 

4150  IF  *.C01-*ABSF*SIGMAoaa4l6'0, 4160, 4161 

4161  IF  %BETA  -  1.04141,4165,4162 
4141  IF  *$IGMAn4166, 4167, 4167 

4162  IF  5SBETA1  -  5.04165,4165,4141 
4160  IF  SSIGMAn  4163,4164,4164 

4163  IF  *B ETA  1—4. 7 1 23895o4 165, 4 1 66, 4 166 

4165  BETA  It  3  .  14  1 59  2  65-  S  I  GMA 
GO  TO  4300 

4166  BETA  It 6.28318530  &  SIGMA 
GO  TO  4300 

4164  IF  XBETAl  -  1.5707963267504167,4165,4165 

4167  BETA  It  SIGMA 

4300  GAMMA  it  AT ANF ** SORT F *  1 . -%ZS  I NF *PHI  Oa «S I  NFSPH I oa G*COSF*PHIOa.COSF*P 

lHIn»C0SF*THETA-THET0noa#*2an/*SINF*PHI0a*SINF*PHIn£C0SFXPHl0n*C0SF 

2*PHIn#C0SF°«THETA-THET0nao 
870  RN  //GAMMA  *  3440.239 

VT1//VT  : 

RA1//  RA  / 

ALT9  It  ALT  ' 

601  CONTR  //CONTR  6  1. 

IF  *CONTR  -l.o4C00,607,604 
604  IF  *CONTR— 9.0602,602,606 
606  CONTR  #0. 

IPAGE  //IPAGE  6  1 
GO  TO  6C1 
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607  PRINT  5017, IPAGE 
602  AAL  //DEL*  57.2958 
A  A2  //BET  A*57. 2958 
A  A3  //THE  T  A*  57 . 2958 
AA4  //PHI  *57.2958 
PRINT  5005,T,X,Y,Z 
AA6  #  SIGMA  *  57.2958 
A A7  #  GAMMA3  *  57.2958 
PRINT  5006, ALT, VX,VY,VZ 
PRINT  5007, VT ,AX,AY ,  AZ 
PRINT  5008, AA2.RN, DEN, GRAY 
IF  SSENSE  LIGHT  In  3000,3001 

3000  PRINT  5028, AAA, AA3,AAl ' 

GO  TO  609 

3001  PRINT  5025, AA4, AA3,AA1,AT 

609  IF  SALT  -  150000. nl940, 1940, 1941 
1540  DELT  //0ELT2 
1941  VX  M  VX  £  AX*DELT 

VY  #  VY  6  AY  *  DELT 
VZ  H  VZ  &  AZ  *  DELT 
X  //  X  £  VX  »  DELT 
Y  I  Y  C  VY  *  OELT 
Z  //  Z  £  VZ  *  DELT 

VT  //  SQRTFSVX**2  6  VY-**2  £  VZ**2n 
AT  //%VT-VTln/DELT 
RA  //  SQRTFSSX**2  £  Y**2  £  Z**2n 
ALT  //RA  -  R 

IF  S ALTn 1202 , 1203, 1203 

1202  A.ABB  ti  A  ABB  £  1. 

1203  IF  SAABB-1.C4000, 2462, 1204 

1204  IF  SAABB-2.c2462.2463, 8CC0 

2463  DEN  # . C 0 2 376 L%% AL T/ %- 5000 . na»* . 0027448-. 002 3769ao 
GO  TO  798 

5028  FORMAT  S1H  ,2X,5HPHI  ,  E 15 . 8 , 2X , 5HTHET A, E 1 5 . 8 , 2 X , SHOE LT A, E 1 5 . 8n 
100  FORMAT  S5E15.8U 

5010  FORMAT  S 1HK , 25X , 42HPHYS ICAL  CONSTANTS  -SAME  FOR  ALL  MISSILES-n 

5011  FORMAT  X1H  ,  25X  ,‘57HGRAV  I  TAT  IONAL  ACCELERATION  AT  SEA  LEVEL  32.174 
1  FT/SEC**2a 

5012  FORMAT  S1H  , 30X , 44HEARTH  ROTATION  2PI/24  X  3600  R ADI ANS/ S ECONDn 

5013  FORMAT  S1H  , 30X , 32HR AD  I  US  OF  EARTH  20,902,890  FEETa 

5014  FORMAT  % 1HK , 5X , 1 2H0 THER  INPUTS^ 

5015  FORMAT  S1H  ,5H0EL  , E 1 5 . 8 , 2 X , 5HBET A  , E 1 5. 8 , 2 X , 5HTHE T A , E 1 5 . 8 , 2X 

1.5HPH1  ,E15.8d 

5016  FORMAT  SIH  , 5HALT  ,  E 15 . 8 , 2X , 5HBALC0 , E 1 5 . 8 , 2X , 5HRA  ,E15.8,2X 

1.5HV0  , E 15 . 8a 

5017  FORMAT  % 1H1 , 79X , 5HP AGE  ,I4n 

8000  PRINT  8001 

8001  FORMAT  % IH1 ,///////// 20 X ,  32H****** **E ND  OF  JOB* ****•**// 1H In 

4999  FORMAT  % 1H1 , 32X , 31HBR0WN  ENGINEERING  COMPANY,  INC . , 4X, 5HPAGE  ,I4n 

5000  FORMAT  S1H  ,30X,35HFREE  FLIGHT  AND  RE-ENTRY  TRAJECTORYn 

5001  FORMAT  S1H  , 35X , 24HPR0GRAMMED  BY  W. 8. WARRENa 

5027  FORMAT  S1H  , 30X , 35HTHE0RE T I  CAL  DERIVATION  -  C.F.OSTNERa 

5002  FORMAT  S1H  ,38X,18H  RESEARCH  STAFF  u 
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5003  FORMAT  *1H  ,3<3X,16H10  NOVEMBER  1962u 

5004  FORMAT  3!  IHK  ,  39X  ,  3 IH  PROGRAM  NO.  SP  9  ENGLISH  UNITS. n 

50050FORMAT  % IHK , 2X , 5HT I  ME  ,E15. 8.2X.5HX  ,  E15 . 8 , 2X , 5HY 

I5H2  ,E15.8a 

5006  FORMAT  % 1H  ,2X,5HALT  , E 1 5 . 8 , 2X , 5HVX  , E l 5 . 8 , 2 X , 5HVY 

15HV2  • E  15 .8a 

5007  FORMAT  *IH  ,2X,5HVT  , E 1 5 . 8 , 2X , 5HAX  ,  E 1 5 . 8 , 2 X , 5H AY 

15HAZ  ,  E 15 . fin 

5008  FORMAT  %\H  ,2X,5HBETA  , £  1 5 . 8 , 2X , 5HR ANGE , E 1 5 . 8 , 2 X , 5HDEN 
15HGRAV  • £15 . fin 


,  6 1 5 . 8 , 2  X , 


,E15.8,2X, 
»E15.8»2X| 
, 615.8, 2X, 


5025  FORMAT  ZUi  ,2X,5HPHI  , E  15. 8 , 2X  , 5HTHE T A , E 1 5 . 8 , 2 X , 5HDEIT A , E 15 . 8,  2X , 
1 5HACCEL , E  15 . fin 
4000  STOP 
END 
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FORTRAN  STATEMENTS 


Metric  Units 


0000000000000  0  0  0  00000000 


CARD  TO  PRINT  80/80 


FREE  FLIGHT  AND  RE-ENTRY 
PROGRAMMED  BY  W.B. WARREN 
FOR  RESEARCH  STAFF 
10  NOVEMBER  1962 


PROGRAM  NO.  SP  9  METRIC  UNITS. 


PHYSICAL  CONSTANTS 

GQ//GRAV  I  TAT  ION  AL  ACCELERATION  AT  SEA  LEVEL  -  9.80665  ME TERS/SEC**2 
U  #  EARTH  ROTATION  -  X2*P lo/^24»3600a  RADIANS/SEC 
R  H  RADIUS  OF  EARTH  -  6,371,230  METERS 
BALCO  #  BALLISTIC  COEFFICIENT  -  CONSTANT. 

MU**2  H  GO *R» * 2  It  9.80665  *  6,371,230**2 

OTHER  INPUT  DATA 

DEL  WANGLE  BETWEEN  HORIZONTAL  AND  MISSILE  AXIS — POSITIVE  UP 
BETA  //AZIMUTH  ANGLE  CLOCKWISE  FROM  NORTH  TO  MISSILE  AXIS 
THETA  //LONGITUDE  OF  MISSILE 
PHI  //LATITUDE  OF  MISSILE 
ALT  //HEIGHT  OF  MISSILE  ABOVE  EARTH  SURFACE 
RA  #  MISSILE  ALTITUDE  FROM  CENTER  OF  EARTH 
VO  H  MISSILE  BURNOUT  VELOCITY  #VT  AT  T#0 
READ  100, DEL, BETA, THETA, PHI, ALT 
R E AC  100, R, VO, GO, U, BALCO- 
READ  ICO , T , OEL  T , DEL  T2 
CON 8  #  0. 

SENSE  L I GhT  1 
CON  1 5  It  0. 

coni  a  c. 

A ACC  #  0. 

CON  TR//0  . 

A  ABB  ii  1. 

BET Al#  BETA  -  .001 
ALT9  #1000. 

RN#C . 

VT  U  VO 
VT1  #V0 

AMU  II  SQRTFSG0*R**2n 

IPAGE  It  1 

BET  AO  #  BETA 

THETO  #  THETA 

PHIC  #  PHI 

RA  U  ALTCR 

AK  U  16 . 0 1 8 A 

OVX  //VO*S  INF^DELn*COSF^PHIa*COSF«THETAa-VO*COSF%DELn*COSFXBETAa*SIN 

lF%PFIn*C0SF$THETAn-V0*C0SF?DELa*SINF3>BETAn*SINF%THETAa 

0 V Y  #VO*S  INF:i!DELn*COSF!SPHIa»SiNFUHE7AQ-VO*COSFSOELn*COSF%BETAa*SlN 

lF%PHIo*S  INFUHETAnCVO*COSF:t;DELn*SlNFSBETAa»COSFUHETAn 

OVZ  #VO*SINF^DELa*SINF35PHIaflvO*COSF%DELn«COSF?BETAn*COSF%PHlP 

X  #RAvCGSF2PHIn*C0SF2THETAa 
Y  #RA*COSF^PH In*SINF%THETAo 
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Z  #  RA*S INFSPH I n 
RA1  #  ALT  G  R 
PRIM  4999,  IPAGE 
PRINT  5000 
PRINT  5001 
PRINT  5027 
PRINT  5002 
PRINT  5003 
PRINT  5004 
PRINT  5010 
PRINT  5011 
PRINT  5012 
PRINT  5013 
PRINT  5014 

PRINT  5015, DEL, BETA, THETA, PHI 
PRINT  5016, ALT,BALCO,RA,VO 
GAMMA  #  0. 

2462  GRAV  #  SG0*R*»2n/RA**2 

IF  SALT  -  ALT9n900,901,901 

900  IF  SCON  1 5  -  l.n3840,901,901 
3840  PRINT  5050 

5050  FORMAT  S 1HK , 3 5X , 6HZ EN I THo 
CON  1 5  #  1. 

901  ALT  n  ALT  *  3.280833 

IF  %ALT- 400000. o4200,4201 ,4201 
4200  IF  SALT  -  40000.nl, 2, 3 

1  D  #  ALT/ 10000. 

R3#  -0 . 3 7728 875 E—  02 
R2#  0 . 5  2352  52  3E-02 

R  1  #  -0.3  1047929EG00 
RO#  0.89444 960EG01 
GO  TO  50 

2  CONTINUE 

3  IF  SALT-80000.n4,5,6 

4  CONTINUE 

5  0  #  SALT-40000. u/IOOOO. 

R3 #  0.23298604E— 05 

R2#  0.21620268E-03 
R1 #-0.479742078600 
RO#  0.754024  19E60I 
GO  TO  50. 

6  IF  SALT  -  160000. n7, 8,9 
70#  SALT-80000. P/10000. 

R3#  0. 2 63 59811 E— 04 

R2  #  0.8  094  3  560E- 02 

Rl#  -0.5 1509227E600 
RO#  0.56272847E601 
GO  TO  50  ' 

8  CONTINUE 

9  IF  SALT-175000. olO, 11,12 

10  CONTINUE 

11  D  #  S  A  L  T-  160000.P/10000. 

R3#-0. 57 066660 E— 04 


CARO  TO  PRINT  80/80 


R2#  0. 30080000E-03 
Ri#- 0.3638 4 173EG00 
RO #  0.20423758EG01 
GO  TO  50 

12  IF  SALT-270000. nl3,  14,  15 

13  D  H  SALT— 175000. a/10000. 

R3#  -0. 84014648 E- 03 

R2#  -0. 16797257E-02 
Rl#  — 0.3293 1657EC00 
RO#  0. 15021266EG01 
GO  TO  50 

14  CONTINUE 

15  IF  SALT-290000. al6,  17, 18 

16  CONTINUE 

17  0  #  SALT-270000. a/10000. 

R3#  0". 

R2 #  0 . 3 6 4 8 0 C 0 9 E- 0 3 
R1#-0.6137624GEG00 
R0#-0.25094S64EG01 
GO  TO  50 

18  IF  SALT-3'50000. ni9, 20,21 

19  0  #  SALT-290000. a/10000. 

R3#  0. 1707C042E-02 

R2  #  — 0.1657S2 18 E—  02 
Rl#  -0. 638987  18EG00 
RO#  — 0.37292327EG01 
GO  TO  50 

20  CONTINUE 

21  D  #  SAL  T- 3500  O'O. a/10000. 

R3#  -0. 30550824 E- 02 

R2#  0.65901945E-01 

Rl#  —  0.687055 5 9EG00 
RO#  — 0.72589563EG01 

50  RO#S EX PFimR3«0GR2a*DGR  la. OGROao/ 100000. 

DEN  #SR0»AK»RA**2a/R*»2 
ALT  #  ALT  /  3.280833 

79  8CAX#2.*VY-*UGU**2*X-:UAMU**2*Xa/RA**3a-*%DEN*GRAV  *VX*VTa/ %2 .  *BALCO 
1  an 

0AY#U**2*Y-2.*VX*U-nAMU**2*Ya/RA**3a-S%DEN*GRAV  *VY*VTn/S2. *BALCO 
Ian 

AZ#S-l.a*imAMU»*2»Za/RA**3aG3SDEN*GRAV  *VZ*VTa/*2. *BALCOoaa 
IF  SC0N8a4CC0, 699, 1569 
699  C0N8  #  C0N8  G  1. 

GO  TO  601 
4201  DEN  #  0. 

ALT  #  ALT  /  3.280833 
GO  TO  758 

1569  IF  SALT  -  457 20 . 096a 1434, 1 434 , 3 104 
1434  IF  SAACCO40C0, 2460,2461 

2460  AACC  #  A ACC  G  1. 

GO  TO  31.04 

2461  CELT  #0ELT2 
3104  T  #  T  G  DELT 


CARO  TO  PRINT  80/80 


422  THET.1  «  THETA 
BETA2  //BETAl 
BETA  1  II  BETA 
SIGKAMSIGMA 
PH  II  U  PHI 

704  PHI  It  ATANFSZ/SSQRTFSRA**2-Z**2oon 
THETA  //ATANFSY/Xo 
OELRA  II  SRA  -  RAla/OELT 

DEL  U  ATANFUDELRAn ASQRTFm**2-DELRA*»2Dnn 
421  DTHET #TH£TA  -THET1 

GAMKA3#  ATANFSSSQRTFSl.-nS  INF%PHI  la*SI  NF2PH I nn6*C0SF*PHI  ln»COSFKP 
lHIn*COSFXDTHETnna##2nn/*SINF*PHIln#SlNF*PHlnCCOSF*PHI  lo*C0SF3!PHla* 
2C0SFSDTHETodo 

SIGMA  H  ATANFS£SINFSDTHETo*C0SFSPHIlaa/SSQRTFS?SINF£GAMMA3oa**2-:?S 
iINFSDTHETa»COSF*PHI loa#*2onn 
IF  SI  .5688  -  SA8SFSSIGMAnon4099, 4160, 4150 
4099  IF  SS  IGMA1Q4102, 4101,4102 

4101  DEtSIG  #  .001 
GO  TO  4120 

4102  DELSIG  §  SIGMA  -SIGMA  1 

4120  IFS6ETA1-BETA204 128, 4130, 4129 

4128  BETA#  BE T A l-S A BSFSDEL S l Gao 
GO  TO  4300 

4130  BETA#  BETA1S.001 
GO  TO  4300 

4129  BETA#  BETA1  GSABSFSOELSIGnn 
GO  TO  4300 

4150  IF  S.C01-XABSF%SIGMAaan4160, 4160, 4161 

4161  IF  SBETA  -  1.04141,4165,4162 
4141  IF  SSIGMAD4166, 4167, 4167 

4162  IF  SBETA1  -  5.04165,4165,4141 
4160  IF  SSIGMAo  4163,4164,4164 

4163  IF  SBETAl-4. 712389504165, 4166, 4166 

4165  BETA  #3.  14  15926 5- SIGMA 
GO  TO  4300 

4166  BETA  #6.28318530  6  SIGMA 
GO  TO  4300 

4164  IF  SBETA1  -  1.5707963267504167,4165,4165  „  . 

4167  BETA  U  SIGMA  '  > 

4300  GAMMA  H  AT ANF USQRTFS  1 . -US  I NF  SPH I  Oo*S I  NF SPH I  antSCOSF SPH  l Oo »COSF XP 
1H  Io*C0SFSTHETA-THET0oon**2no/*SINF%PHI0a*SINFSPHIn£C0SF!SPH  IOa*COSF 
2SPHIn*C0SFSTHETA--THET0onn 
870  RN  //GAMMA  *  3440.239 
VT1//VT 
RA1#  RA 
ALT9  HALT 

601  CONTR  //CONTR  i  1. 

IF  SCONTR  -1.04000,607,604 
604  IF  SC0NTR-9. o602, 602, 606 

606  CONTR  #0. 

IPAGE  U IPAGE  €  1 
GO  TO  601 

607  PRINT  5017, IPAGE 
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CARD  TO  PRINT  80/80 


602  AAl  It  DEL*  57.2958 
AA2  //BETA*57. 2958 
A A3  «THETA*57.2958 
AAA  It  PHI  *57.2958 
PRINT  5005»T»X»Y»Z 
A A6  «  SIGHA  *  57.2958 
A A 7  #  GAMMA 3  *  57.2958 
PRINT  5006,  ALT,  VX,VY,VZ. 

PRINT  5007, VT,AX,AY ,AZ 
PRINT  5008  , AA2,RN, DEN, GRAV 
IF  SSENSE  LIGHT  In  3000,3001 

3000  PRINT  5028, AAA, AA3, AAL 
GO  TO  609 

3001  PRINT  5025, AAA, AA3, AAl, AT 

609  IF  SALT  -  A5720.096O19A0, 19A0, 19A1 
1 9 A 0  DELT  It 0ELT2 
1 9 A 1  VX  it  VX  0  AX*DELT 

VY  It  VY  8  AY  *  DELT 

VZ  A  VZ  £  AZ  *  DELT 

H  X£  VX  »  DELT 

Y  #  Y  6  VY  *  DELT 

Z  #  Z  S  VZ  *  DELT 

VT  V/  $QRTFSVX**2  £  VY**2  8  VZ**2n 

AT  f/SVT- VTln/DELT 

RA  /f  SQRTF%X**2  8  Y**2  £  Z**2a 


ALT  ilR A  -  R 

IF  SALTD1202, 1203, 1203 

1202  AABB  It  A  ABB  8  1. 

1203  IF  SAA8B-l.nAC00,2A62,120A 
1 20 A  IF  SAABB-2.n2A62,2A63,8CC0 

2A63  DEN  it . 00 23 768SS AL T/ %- 5000. on* % . 0027 A A8-. 002 3769an 
GO  TO  798 

5028  FORMAT  S1H  ,2X,5HPHI  ,  E  1 5 . 8 , 2X  ,  5HTHE  TA ,  E 15 . 8 , 2  X  ,  5HDELT  A ,  E'l5 . 8n 
ICO  FORMAT  S5E15.8n 

5010  FORMAT  % 1HK , 25X , A2HPHYS ICAL  CONSTANTS  -SAME  FOR  ALL  MISSILES-a 

5011  FORMAT  S1H  ,  25X , 62HGRAV I  TAT  I ONA L  ACCELERATION  AT  SEA  LEVEL  9.8066 
15  METERS/SEC**2u 

5012  FORMAT  S1H  , 30X , AAH EARTH  ROTATION  2PI/2A  X  3600  RADI ANS/SECONDo 

5013  FORMAT  SlH  , 30X ,  33HRA0I US  OF  EARTH  6,371,230  METER  Sn 

501  a  format  sihk, 5x, 12hother  inputsd 

5015  F  u  ,<  Mi  A  T  S1H  ,  5HDEL  ,  E  1 5 . 8 , 2  X ,  5HBE  T  A  ,  E 1 5 . 8 , 2X  ,  5HTHET  A ,  E  1 5 . 8 , 2X 


1.5HPHI 
5016  FORMAT 
1 , 5HV0 


, E15. 0n 

% 1H  , 5HA  LT  ,E15.0,2X,5HBALCO,E15.8,2X,5HRA 
,El5.8n 


,  E  15i>8 , 2X 


5017 

FORMAT 

S1H1 

8C00 

PRINT 

8001 

8001 

FORMAT 

S1H1 

AS99 

FORMAT 

%  1H1 

5000 

FORMAT 

%  1H 

SCO  1 

FORMAT 

%  IH 

5027 

FORMAT 

S1H 

5002 

FORMAT 

*1H 

5003 

FORMAT 

X  IH 
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CARO  TCJ  PRINT  80/80 


5004  FORMAT  % 1HK , 39X , 3 LHPROGRAM  NO.  SP  9  METRIC  UNITS. n 
,  50050FORM AT  %  IHK  »  2X • 5HT I  ME  , E15 . 8 , 2X  , 5HX  , E15 .8,2X ,5HY  ,E15.8,2X, 

15HZ  ,E15.8n 

5006  FORMAT  HUH  .2X.5HALT  , El  5 . 8 , 2X  , 5HVX  ,  E15 . 8, 2X , 5HVY  ,EL5.8,2X, 

15HVZ  ,E15.8a 

5007  FORMAT  $1H  ,2X,5HVT  , E l 5 . 8 , 2X , 5HAX  , E15 . 8 , 2 X , 5HAY  ,E15.8,2X, 

15HAZ  ,EI5.8o 

5008  FORMAT  tlH  ,2X,5HBETA  , E l 5. 8 , 2X .5HRANGE , E 15 . 8 ,2X , 5H0EN  ,E15.8,2X, 
15HGRAV  ,E15.8a 

5025  FORMAT  *1H  ,2X,5HPHI  , El  5. 8 , 2X  .5HTHETA , E15 . 8 ,2X , 5HDELT A, E 15.8, 2X , 
1 5HACCEL , E 15 . 8n 
j  4000  STOP 

I  END 

II 


[ 


ri 

LJ 


."SPTM 

udiJ 


rmn 

Jfc 

1 

1 

I 
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